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A first-principle study of tfie collective oscillation spectrum of a strongly correlated one-component 
plasma in a strong magnetic field is presented. The spectrum consists of six fundamental modes 
which are found to be in good agreement with results from the Quasi-Localized Charge approxima- 
tion. At high frequencies additional modes are observed that include Bernstein-type oscillations and 
their transverse counterparts which are of importance for the high-frequency optical and transport 
properties of these plasmas. 
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Strongly coupled charged particle systems in which 
the interaction energy exceeds the kinetic energy are of 
rapidly growing interest in many fields reaching from con- 
densed matter to ultracold gases, trapped ions, nonideal 
plasmas up to the quark-gluon plasma, e.g. |1J. Quite of- 
ten these systems are subject to a strong magnetic field - 
examples are white dwarf stars, neutron star crusts OH], 
magnetized target fusion scenarios [4[ |5j or quantum 
Hall systems [6 . Magnetic field effects are also com- 
ing into the focus of laboratory experiments with dusty 
plasmas [1], ultracold neutral plasmas [7] and trapped 
ions [8]. Apart from basic transport properties such as 
diffusion [9], the collective excitation spectra are of fun- 
damental relevance e.g. for the response to optical ex- 
citation. In two-dimensional (2D) plasmas, the prop- 
erties of these wave spectra are firmly established ana- 
lytically via the Quasi-Localized Charge Approximation 
(QLCA) of Kalman and Golden p!QHT2] and computer 
simulations [13 . Strong correlations lead to the exis- 
tence of a shear-type mode |T5] alongside the famil- 
iar compressional modes (plasmons) known from weakly 
coupled situations. These shear modes could be verified 
in 2D dusty plasma experiments [I6l |T7] . In the pres- 
ence of a strong magnetic field, these transform into two 
modes commonly known as magnetoplasmon and -shear 
or upper and lower hybrid mode, e.g. [18 . 

Surprisingly, the corresponding behavior of the much 
more common strongly coupled three-dimensional plas- 
mas has remained unexplored, except for a preliminary 
study in Ref. [19 . The present work is, therefore, de- 
voted to filling this gap by a theoretical approach based 
on the QLCA and complementary first-principle molec- 
ular dynamics (MD) simulations. We obtain the com- 
plete set of fundamental linear oscillations and explore 
their wave-vector dispersion and polarization. Further- 
more we present evidence for the existence of nonlin- 
ear modes that are generalizations of Bernstein modes 
of high-temperature plasmas [20l |2T] . 

Model and simulation method. We consider a clas- 
sical one-component plasma (with a neutralizing homo- 



geneous background, OCP) subject to an external mag- 
netic field B = Bbz' The exact equations of motion for 
N particles are 

= F^/m + cJcVi X e^, i = l...N, (1) 

where = and Fi is the Coulomb force due to all 
particles j i. The thermodynamic equilibrium state 
of the system is characterized by the coupling parame- 
ter T = X {AireoakBT)-^, where a = [3/ (4n7r)]^/^ is 
the Wigner-Seitz-radius with the number density n, and 
q denotes the particle charge. The strength of the mag- 
netic field B is given by /3 = ujc/oop ex 5, i.e., the ratio 
of the cyclotron frequency, ujc = qB/m^ and the plasma 
frequency, cj^ = nq^ / {som). The spectra of the longitu- 
dinal and transverse fluctuations, Lik^uo) and T{k^uj)^ of 
the OCP follow from the Fourier component of the gen- 
eral current density operator in standard manner [22l [23] 

N 

j(/c,t) = ^v,(t)exp[zk.r,(t)], (2) 
via a temporal Fourier transform J't 

^}}^J\MKk,t)}f. (3) 

The current can be decomposed into a part parallel and 
a part perpendicular to k, j = jy + jx, so application of 
(|3| to j|| yields the longitudinal spectrum L{k^uj) and jx 
tne transverse spectrum T{k^uj). Collective oscillations 
appear as peaks in these spectra, cf. Fig. [l] 

Zero magnetic field. To set the stage, consider first 
the unmagnetized system. In Fig. [l] we plot the longitu- 
dinal and transverse and longitudinal fluctuation spectra 
L{k^uj) and T{k^uj) together with the well-known QLCA 
results for the collective modes (solid black lines p!4| l24 ] ). 
There are two remarkable deviations from the familiar 
spectrum of a weakly coupled plasma, e.g. [20l[23: first, 
the plasmon dispersion does not show a monotonic in- 
crease as cj^(/c) = ujp{l + 3/c^r|)) [rjj is the Debye radius] 
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FIG. 1. (color) Longitudinal spectrum L(k,uj) and trans- 
verse spectrum T{k^uj) of an unmagnetized OOP at F = 100. 
The black lines are the QLCA dispersion relations (Plasmon 
and ordinary shear mode) for F = 100 (solid line) and F = 20 
(dashed line). 



but decays and exhibits oscillations. Second, there ex- 
ists an additional shear mode in the transverse spectrum. 
Fig. [l] also shows that the precise value of F is of minor 
importance, increase of the coupling tends to slightly in- 
crease the amplitude of the oscillations. Overall, Fig. [l] 
indicates that QLCA describes the oscillation spectrum 
of a strongly coupled OCP rather well [26 . 

Theory for a magnetized OCP. We now turn to a 
theoretical description of the wave spectrum in a magne- 
tized strongly correlated OCP using the QLCA approach, 
e.g. The equation of motion for the collec- 

tive coordinates ri\^a{^) in the presence of a magnetic 
field read 



^2 + i(^(^cCra/3 ^ka(^) = 0, (4) 
where a 




and Do^p is the QLCA dynamical matrix which is 
a functional of the pair distribution function g{r). The 
non-trivial solutions to Eq. Q represent the modes of 
the system which turn out be rather complicated. There- 
fore, here we concentrate on modes with k || B or k ± B. 
Their main properties are summarized in Tab. [T| and the 
dispersions are depicted for three different magnetic field 
strengths at strong coupling (F = 100) in Fig. [5] Obvi- 
ously, there exist five possible (pairwise parallel or per- 
pendicular) orientations of the vectors B,k and r - two 
(three) of them corresponding to k || B (k ± B), cf. 



Tab. [l| At the same time, the QLCA equation Q yields 
six fundamental solutions - three for each orientation of 
the wave vector - which we discuss in the following. 

Wave vector parallel to B: The three solutions of 
this type (black lines in Fig. [2]) are the longitudinal plas- 
mon (P) which is unaffected by the magnetic field, since 
the particles oscillate parallel to B, and two shear-mode 
solutions (upper and lower shear, US, LS). The latter 
arise from the shear mode of an unmagnetized plasma, 
cf. Fig. [l] whose degeneracy is lifted by the magnetic 
field. Their frequency difference equals Uc (cf. Tab. |l| 
and thus grows linearly with B. 

Wave vector perpendicular to B: The three solu- 
tions (red lines in Fig. |2| are the ordinary shear mode 
(OS) and the upper and lower hybrid modes (UH and 
LH) . The ordinary shear mode is ^-independent since the 
particle oscillation occurs along B. The long- wavelength 
limit of the UH (LH) mode is uj = {uo^ + a;2)V2 = q). 
These two modes have peculiar polarization properties. 
While all four other modes are characterized by a fixed 
orientation of the displacement vector relative to k - 
there are one longitudinal (P) and three transverse os- 
cillations (US, LS, OS), the UH and LH modes ex- 
hibit particle oscillations that rotate with time ellipti- 
cally in the plane perpendicular to B. To clarify these 
properties, the bottom part of Fig. |2] shows the mod- 
ified eccentricity e = — /a^^ where the semi- 
major axis a = max{|?7a;|, 1%!} and the semi-minor axis 
b = min {\r]x\AVy\} 7 equals 1 if r^^^ > r]y and — 1 
otherwise. Here, fj is the eigenvector of Eq. Q for the 
solutions UH and LH, respectively. £ = +1 (—1) for lon- 
gitudinal (transverse) waves. Moreover, the shape of this 
ellipse is /c-dependent, cf. Fig. [2] [29]. The UH mode 
starts longitudinal, at long wavelengths, and transforms 
into a predominantly transverse wave, at the first cross- 
ing of the P and OS modes. At each subsequent cross- 
ing point, the dominant character of the wave changes 
again. The LH mode exactly mirrors this behaviour. For 
/c ^ 00, both modes are in-plane circular. The limit 
B ^ of these modes is particularly interesting: while 
all 6 modes converge to the two modes shown in Fig. [l] 
the LH and UH modes only do so in a pointwise manner. 
The LH, at all /c, remains below the UH, i.e. we observe 
an avoided crossing of the two at the intersection points 
of the plasmon and shear mode (see the case /3 = 0.1 in 
Fig.|2|. 

Simulation results. Let us now critically test the 
theoretical predictions against the computer experiment. 
To this end, Eq. ([T]) is solved by a microcanonical MD 
method for TV = 8000 particles, incorporating an arbi- 
trarily strong static homogeneous magnetic field into the 
second-order Velocity Verlet algorithm ^ . Prior to data 
acquisition, the system is equilibrated by an isokinetic 
thermostat. The above mentioned five possible orienta- 
tions of B,k and r manifest themselves in five different 
currents of the type ([2| leading to five fluctuation spec- 
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Polarization 


Dispersion uJi{k) 


Asymptotic uo^ 


Spectrum 


Plasmon (P) 
Upper Shear (US) 
Lower Shear (LS) 


II 
II 

II 

II 


II 
II 




longit . 
transv. 
transv. 


circular 
circular 


y^p -r JyL[K) 

^[y/uJ^,+4DT{k)-UJc] 


(jOE 

L 




J 1 


Ordinary Shear (OS) 
Upper Hybrid (UH) 
Lower Hybrid (LH) 




II 


t-dep. 
t-dep. 


transv. 
hybrid 
hybrid 


linear 
in-plane elliptical 
in-plane elliptical 


VDrik) 
^\^Ql + DT{k) + Q^^ 
^[Qi + DT(^)-n^]' 


ooe 

\ y^ul + 4cj| + 
\ [7cj2 + 4^1 _ 







TABLE L The six principal collective modes in a magnetized three-dimensional OOP in QLCA. Shown are the relative orienta- 
tions of the wave vector /c, magnetic field and particle displacement r, the wave type, the polarization, the dispersion relation 
0Ji{k), its asymptotic, cj°°, for large k and the current spectrum in which they appear, uje = ujp/y/S is the Einstein frequency, 

Dl/t denote the longitudinal and transverse parts of Da/3, — 
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FIG. 2. (color online) The modes of a 3D magnetized OOP 
(/^ = 0, r = 100) in QLCA, Plasmon (P), Ordinary shear 
(OS), upper hybrid (UH), lower hybrid (LH), lower shear 
(LS) and upper shear (US), for three magnetic field strengths. 
Modes with k || B are shown in black (dark), and those with 
k ± B in red (light). Bottom: /c-dependence of the eccentric- 
ity e of the UH (black) and LH (red) for (3 — 0.1 (solid) and 
P = 0.5 (dashed). 

tra: there are two longitudinal ones, and [the 
superscript denotes the mutual orientation of B and r] 
and three transverse ones, T", T,!' and t\_ [the subscript 



specifies the angle between B and k in case of an ambigu- 
ity]. The correspondence of the five fluctuation spectra 
to the six collective oscillations can be seen in Tab. [H 

Fig. [3] shows the wave spectra of a magnetized OCP at 
F = 100 and /3 = 1.0. To capture the two hybrid modes, 
the two spectra L^(/c,cj) and T^(k^uj) have been added 
together. Six frequency peaks can be identified which 
overall show good agreement with the QLCA modes. 
This conclusion is representative and holds also for other 
values of F and j3. At the same time, there exist three 
noticeable deviations of the QLCA from the simulations. 
The first is observed at low frequencies (corresponding 
to long times) due to the breakdown of the QLCA as- 
sumption of a frozen potential landscape [12]. Second, 
the low-frequency shear modes vanish at small /c, which 
is not reproduced in QLCA ^T] [32] . Finally, it has been 
observed that QLCA misses the high-frequency part of 
the spectrum in a 2D magnetized OCP where Bernstein- 
type modes were found in the simulations [33l |34]. It 
is, therefore, tempting to analyze more in detail the fre- 
quency region that is not captured by Fig. [3| 

In the left column of Fig. [4] we present the high- 
frequency part of the spectra corresponding to B ± k 
at three values oik. In addition to the two known modes 
(UH, LH), the spectrum + exhibits a number 
of further peaks that are equally spaced. While these 
peaks are strongly damped for /3 = 1, they become very 
sharp at a five- fold stronger magnetic field, cf. center 
column of Fig. [4| These modes in the spectrum of waves 
propagating perpendicular to the magnetic field clearly 
resemble Bernstein modes of high-temperature plasmas, 
e.g. [20l El]. In a 2D OCP it was found recently that 
"dressed" (modified by correlations) Bernstein modes 
reappear at strong coupling [33, 34 . The analysis of 
the present 3D spectra shows that these modes appear 
at multiples of the high-Zc limit of the upper shear and 
upper hybrid frequency (see the asymptotics in Tab. |l|. 
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FIG. 3. (color) The complete set of oscillation spectra of 
a magnetized OOP at F = 100 and /3 = 1.0 for frequencies 
uj < 2ujp. The solid lines correspond to the six QLCA modes 
that are labeled in the figure, cf. Tab.[l] The orientation of the 
magnetic field (blue), wave vector (orange) and the particle 
oscillations (red) are indicated by the arrows. 



^UH = ^US = ^oc = Ui'^c + 4cj|)l/2 + CJc], with the 

Einstein frequency Ue = ujp/\^. 

Apart from this similarity to a 2D OCP, the high- 
frequency spectrum in 3D is much richer, and the har- 
monics nujoo show up also in other fluctuations. Indeed, 
as can be seen in the lower left part of Fig. |4] the har- 
monics emerge also in the spectrum T" around the same 
frequencies (solid grey lines) and more clearly at stronger 
fields, cf. bottom center figure. A closer inspection 
shows that there these peaks are visible only at small ka 
whereas, at higher /ca, instead, a double peak structure 
emerges at nco'^icj^;. This is accompanied by a broaden- 
ing of the higher harmonics peaks in L^(cj) +T^(cc;) that 
develop "shoulders" around nuj oo :^oje- A similar trend 
is observed in the two B || k excitation spectra, cf. right- 



most column of Fig|4] Here, the US mode is clearly seen 
in T||^(cc;), as a peak around u/up ^ 5. Additionally, the 
second harmonic of the US is visible around uj/up ^ 10. 
These two peaks also appear in (cj), alongside the low- 
frequency plasmon peak ai uo/uOp ~ 1/a/3, cf. Tab.jlj The 
peaks at cj/cjp ~ 5, 10 are accompanied by the same dou- 
ble peak structure seen in the cross-field wave spectra. 

All these high-frequency features can be understood 
within a common simple physical picture: scattering of 
multiple short wavelength oscillations where uji -\- UJ2 = 
uji2 and ki +k2 = ki2, a mechanism well-known in high- 
temperature plasmas [20l [25] and nonlinear optics. In 
fact, the n-th harmonic of the US (UH) is explained by a 
cascade of nonlinear inelastic interactions c^oo ^ ^cjoo 
. . .nuJoo' On the other hand, the pair of side peaks at 
UJ oo ± in the middle column of Fig. |4] arise from the 
scattering of an UH (frequency cjoo) with an OS wave 
(frequency Ue)^ and the analogous peaks in the right 
column of Fig. [4] are due to a scattering of an US mode 
and a plasmon (frequencies c^oo and oje)- Finally, the 
appearance of the peaks at uj/ujp ~ 5,10,15 in T" [at 
uj/ujp ^ 5,10 in LII] are due to elastic scattering of an 
UH [US] mode (or a corresponding higher harmonic) with 
that frequency on the particles during which the particle 
oscillation direction is turned from in-plane (± B) to out- 
of plane (|| B) which requires a finite shear elasticity - a 
remarkable feature of strongly correlated plasmas [35 . It 
is particularly striking that this mechanism also excites 
novel longitudinal modes in the plasmon spectrum 
which have a 5-dependent frequency, in contrast to the 
fundamental plasmon. 

In summary, we have presented the complete collec- 
tive oscillation spectrum of a strongly correlated one- 
component plasma in a magnetic field. The OCP sup- 
ports six base modes - three for each k - which are well 
described by the QLCA, and their characteristics have 
been given in Tab. [Tj In addition to these modes there 
exist additional oscillations arising from elastic and in- 
elastic wave scattering. They are fuelled by the familiar 
Bernstein scenario of harmonics generation in a magnetic 
field. However, in contrast to high-temperature plasmas 
the base frequency differs from Uc - it is modified by cor- 
relations. Furthermore, correlations give rise to a novel 
class of modes: transverse Bernstein waves which are 
harmonics of the upper shear oscillation. The present 
results are representative for a broad range of coupling 
parameters and for Coulomb and Yukawa OCP as well. 
They are, therefore, expected to be of direct relevance 
for the high-frequency transport and optical response of 
a broad spectrum of strongly coupled magnetized plas- 
mas - from white dwarf and neutron stars to trapped 
ions and dusty plasmas. 

This work is supported by the DFG via SFB-TR 24 
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gram, and a grant for CPU time at the North-German 
Supercomputing Alliance HLRN. 
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FIG. 4. (color online) Fluctuation spectra of an OOP extended to high frequencies, at two values of the magnetic field and 
three wave numbers, as indicated in the figure. Vertical lines indicate the harmonics ncjoo and the frequencies nuoo =t cj^;. 
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